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Superconductivity in the quasi-one-dimensional material Lio.gMo6 0i7 is analyzed based on a mul- 
tiorbital extended Hubbard model. We found strong charge fluctuations at two different momenta 
Qi and Q2 giving rise to two different charge ordered phases. Evaluating the superconducting 
vertex, we found superconductivity near strong charge fluctuations at Qi. The order parameter 
has the p-wave symmetry with nodes on the Fermi surface. The metallic state displays a char¬ 
acteristic charge collective mode Qi due to nesting and for on-site Hubbard repulsion sufficiently 
large, a charge critical mode Q 2 driven by Coulomb repulsion, which softens at the proximity to the 
transition. The results are quite robust for different coupling parametrizations. A phase diagram 
discussing the relevance of the model to the physics of the material is proposed. 

PACS numbers: 71.10.Hf, 71.10.Fd, 74.40.Kb, 74.70.-b 


I. INTRODUCTION 

Low-temperature physics of correlated materials is of¬ 
ten characterized by the competition between ordered 
phases and unconventional superconductivity. Tipically, 
a static mean-field description, implying negligible fluc¬ 
tuations beyond the limits of the ordered phase, is not 
valid in these systems. Nearly all dynamical probes 
show strong order parameter fluctuations, not only in 
the neighboring superconducting phases, which suggests 
a natural mechanism of pairing, but also in the strange 
metal, present at higher temperatures. Lithium Pur¬ 
ple Bronze (LiPB), adds the ingredient of quasi-one- 
dimensionality to the problem and suggests the possi¬ 
bility that charge and spin fluctuations alone, without 
the existence of real order, might be responsible of su¬ 
perconductivity and anomalies of the normal phase. 

The metallic phase of LiPB, with chemical formula 
Lio.gMo60i7, has been characterized as a robust Lut- 
tinger Liquid (LL) in a series of Angle Resolved Pho¬ 
toemission Spectroscopy (ARPES) experiments rang¬ 
ing different temperature regimes, sample growth tech¬ 
niques, photon energies, and data analysis procedures 1 9 . 
STM spectroscopy shows 10 11 LL single-particle density 
of states and thermal and electric transport measure¬ 
ments are in complete disagreement with Widemann- 
Franz law 12 . When temperature is decreased, an up¬ 
turn of the resistivity occurs at T m ~ 20 K 12 15 and the 
material becomes superconducting at lower temperatures 
around T c ~ 1 K ! 14 * 10 [ 

Unlike other low-dimensional bronzes, the resistivity 
upturn of LiPB 17 is not associated with a lattice dis¬ 
tortion (See Table 1 in Ref. IT8T) . Neither thermal 
expansion 19 nor neutron scattering experiment^ 0 have 
identified a phase transition at T m suggesting the idea of 
a soft crossover of electronic nature. No gap has been 
clearly observed in the spectroscopies but optical con¬ 


ductivity measurement^ suggest the presence of a weak 
pseudogap. Recently, thermopower 22 and NMR 23 exper¬ 
iments have confirmed different aspects of the quasi-one- 
dimensionality of this material but the nature of the up¬ 
turn remains a mistery. 

The most recent study of superconducting properties 16 
confirms quantitatively that the large anisotropies ob¬ 
served in the upper critical field agree with those 
expected from the electrical resistivity in the metallic 
phase. The coherence lengths perpendicular to the 
chains are larger than interchain distances and H c 2 
increases monotonically with decreasing temperature to 
values 5 times larger than the estimated paramagnetic 
pair-breaking field. Neither spin-orbit scattering nor 
strong-coupling superconductivity seem to explain 
this behavior suggesting the possibility of spin triplet 
superconductivity. A quantitative comparison with 
experiment^ shows that superconductivity can be 
destroyed through orbital effects at fields higher than 
the Clogston paramagnetic limit provided that the 
superconducting pairs are in the triplet state. 

In the last years there has been a very important the¬ 
oretical effortPIHSZ] to reduce the complexity of the unit 
cell to microscopic Hamiltonians reproducing different as¬ 
pects of this phenomenology. In this article, we present a 
microscopic theory for the unconventional superconduct¬ 
ing properties observed in Lio.gMo60i7. Based on a min¬ 
imal extended Hubbard model introduced in Ref. I24|25[ 
we show that Lio.gMo 60 i 7 superconducts in the triplet 
channel when charge and spin fluctuations are enhanced, 
which may be also related with the upturn in resistivity 
at T m ^. Using the random phase approximation (RPA), 
we identify the CDW pattern characterized by two or¬ 
dering wave vectors, Qi and Q 2 . In the proximity of 
those phases we evaluate and analize the superconducting 
vertex finding dominant p-wave triplet superconductivity 
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FIG. 2: (Color online) (a) Fermi surface with two bands, Qi a 
nesting vector, and Q 2 referred in the text.(b)Real part of the 
bare susceptibility in momentum space for uj = 0 and q a = 0. 
Notice the maximum reveals the warping of the Fermi surface 
at the nesting vector. 


FIG. 1: (Color online) Schematic crystal structure of 
Lio.gMo 60 i 7 projected onto the b-c plane showing only the 
partially filled Mo atoms forming the zig-zag ladders relevant 
to the low energy electronic properties. Our choice of unit cell 
is highlighted and the orbitals numerated (solid line) accord¬ 
ing to the text, the hoppings (dotted line) and the Coulomb 
interactions (dashed line) are also represented. 

with nodes on the Fermi surface. Within our methodol¬ 
ogy we find results compatible with the one presented in 
a very recent preprint 7 

II. MICROSCOPIC MODEL. 

The electronic structure close to the Ep and the quasi- 
one-dimensionality of the system derives from two par¬ 
allel zig-zag Mo-0 chains per unit celP^ Fig. [l] Tight 
binding® and dftEH band structure calculations agree 
that the Mo-0 orbitals of the chain give rise to four 
bands and two of them cross the Fermi level. ARPES 
confirms the quasi-one-dimensionality of the Fermi sur¬ 
face. A Slater-Koster tight binding parametrization of 
the system has been propoposed in Ref. [24] and the role 
of long-range Coulomb couplings in the anomalies of the 
metallic phase has been also studied 27 . Here, we con¬ 
sider a strongly correlated model, which can capture the 
essential physics of Lio.gMoeO^ 2 ^ consisting on an ex¬ 
tended Hubbard lattice with 4 Mo-atoms per unit cell, 
which reads: 


The diagonalized Hamiltonian: Hq = 

Xlk/icr e M(^ c )^k/icr^ k Mc r ’ leads to four bands denoted 
by /i, the two lowest ones cross the E ^ 24 1 27 I 3 L 1 The 
Fermi surface, close to one quarter-filling, n = 0.225, is 
shown in Fig. [ 2 ] (a). 

The Coulomb interaction terms in the Hamiltonian in¬ 
cludes on-site Hubbard interaction (£/), intra-ladder in¬ 
teraction with the following non-zero matrix elements: 
Vi2 = V32 = V\\ and Vi 2 = V34 = V± and inter-ladder 
W interactions: W 13 = W and W 12 = W34 = Wj_, as is 
shown in Fig [l] (dashed cell). 

Z,z,a l,i,a,j,/3 

+ £ 

( 2 ) 

The interacting Hamiltonian only includes density- 
density Coulomb interaction contributions. Within this 
work, we have consider several combinations of param¬ 
eters, all of them leading to essentially the same results 
presented here where we reduce the parameter space to 
two variables (U and V). We take the Coulomb in¬ 
teraction among different sites with l/|r| dependence, 
where |r| is the distance among orbitals. Therefore, we 
parametrize the interactions by weighting the V’s with 
the interatomic distances: V = Vj|r|| = V±r± = Wrw = 

W±r w ±- 


H = H 0 + H V , (1) 

where Hq is the non-interacting tight-binding Hamilto¬ 
nian. The one-electron Hamiltonian can be expressed in 
terms of Bloch waves with the following non-zero ma¬ 
trix elements, the intra-ladder: t\ 2 (k) = £43 (k) = t± = 
—0.024 eV,and £14 = £ 23 (k)£ = 0.5 eV and the hoppings 
among chains: £13 (k) —t' — 0.036 eV, as is shown in Fig 
[l] (dotted cell). 


III. MULTIORBITAL RPA APPROACH 

In this section we explain the multi-orbital random 
phase approximation (RPA) approach for this model, we 
will study spin and charge ordering based on spin and 
charge susceptibility respectively, and the superconduct¬ 
ing vertex based on projections of different order param¬ 
eters. 
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A. Spin susceptibility 

The RPA spin susceptibility readd^ 

(xs) Q ,/3(q) = (xoW(q) 

+ E(^)«'/3'(q)(^) a ' /3 '(Xo)a /3 (q) (3) 

a' /3' 


where the indices <a, /? refer to the four Mo d xy orbitals 
present in the unit cell. This is the more general case 
for density-density interactions. In our case the spin in¬ 
teraction is a diagonal matrix ( U s ) a p = U5 a ,p, momen¬ 
tum independent. The non-interacting susceptibility, xo, 
reads: 


(Xo)a,/9(q,*w) = 




N ^ 

k ,H,V 


a g( k ) a ^*( k ) fl ^( k + q) a r( k + q) 

w + e„{ k + q) - e M (k) 


[/(e„(k + q)) - /(e M (k))], 


(4) 


where N is the number of lattice sites, and z/, (i are band 2V(q) a ^ where V(q) is the Fourier transform of Vi a jp 

indices. The matrix elements a^(k) = (a\fik) are the and Wi a jp interactions in real space, 
coefficients of the eigenvectors diagonalizing H 0 . 


B. Charge susceptibility 

The RPA charge susceptibility readd 32 : 

(Xc)a,/?(q) = (Xo)a,/3(q) 

- E(^)^'(q)(^c) a ^(q)(Xo)^(q) (5) 

a' (3' 

Where U c ( q) is the Coulomb matrix appearing in Eq. 
(J2| expressed in momentum space: (I/ C ) a ^(q) = + 


C. Superconducting Vertex 

Assuming that the pairing interaction arises from the 
exchange of spin and charge fluctuations, we can calcu¬ 
late the pairing vertex using the RPA, (For a detailed de¬ 
scription of the method, see for instance 32 ). The strength 
of the interaction is weighted by and making use of 
the Kramers-Kronig relation we only need the z ero fre¬ 
quency vertex,^!. For the multiorbital cas d 35 ^ 4 ! , singlet 
and triplet pairing vertex at zero frequency are given by: 


r a3 glet (k, k') =(u+ ^U sXs (k - k')U s + V(k - k') - ^-U c { k - k')x c (k - k')U c (k - k') 


a(3 


r lT let (k, k') = ( -^f/,x s ( k - k')U s + V(k - k') - ^U c (k - k')Xc(k - k’)U c (k - k ') 


(6) 

(7) 


a/3 


We transform the vertex in real space a/3 into momentum an outcoming momentum of (k 7 ,—k 7 ). We take the sym- 
space fiu with the band structure eigenvalues (k) . The metric and antisymmetric parts for singlet and triplet 
Cooper pairs have an incoming momentum of (k,—k) and channels respectively. 


rt nglet (k,k') = ^a“*(-k)a“*( k )Real [Cf et (k,k')l a£(k>£(-k') + (k' o -k') 

a/3 


a/3 


(8) 


r«P let (k,k') = ^>“*(-k)a«*(k)Real 

a/3 


rlt let (k,k') a^(k')a^(—k') - (k' ^ -k') 


(9) 
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FIG. 3: (Color online) (a) Phase Diagram U — V. Extended 
region of p y superconductivity with nodes in the Fermi surface 
(inset) close to the CO region. In the inset we show d x 2_ y 2 
and p y wave functions. 



Eig Re% c 11=0.4 V=0.7 



FIG. 4: (Color online) Left: momentum space distribution of 
the interaction, we show the sum of all components. Right: 
The same figure as Fig. [ 5 ] (bottom right), in those momentum 
relevant to the pairing vertex. 


We solve the gap equation by projecting out and / 
waves P): 


A^y - 


|v F (k^)| IfS Ft(kD|#7(Kx)rJL( 




( 10 ) 

where 7 numerates the different waves projected (s,p,d or 
/) and P depends on the 7 symmetry. P could be singlet 
or triplet. The gap equation has a solution when A 7 is 
1 . We increase the interaction parameters until the dom¬ 
inant wave solves the equation, for stronger interactions 
the gap is already opened in that channel. 


IV. PHASE DIAGRAM 

Using the parametraization described in section [TT] we 
can study the complete parameter space, reduced to two 
variables U and V. The RPA spin susceptibility (Eq. 
[3| breaks at U = 0.47 indicating a Spin Density Wave 
(SDW) phase. 


FIG. 5: (Color online) Top. Imaginary part of the larger 
eigenvalue of the charge susceptibility near the critical value 
(See Fig. |3j for t/ = 0.2 and (U = 0.3) in left and (right) 
panel. Bottom. Real part of the larger eigenvalue in momen¬ 
tum space and zero frequency close to the critical value for 
U — 0.2 and (U = 0.3) in left and (right) panel. 



FIG. 6: (Color online) Gap squared of the Q 2 critical mode 
scaled with the critical interaction V = V c . We observe a | 
exponent near the critical value. 


The RPA charge susceptibility (Eq. [5| diverges for 
different momenta for different U on-site Hubbard 
interaction, leading to different charge order regions 
in the phase diagram (See Fig. [3|. The charge order 
susceptibility divergence consists in an interplay between 
the bare susceptibility strongly peaked at « 7 r (Fig. 
§b)) and the charge interaction U c ( q). The analysis 
involves 4x4 terms but in essence can be understood 
with the sum of the 16 contributions. We observe 
that J7 c (q)(Fig. [4] c) is minimum at the (27r, 2 tt) edge 
of the Brillouin Zone, notice that the periodicity is 
not required since we are dealing with the sum of the 
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elements of a matrix. For U = 0, red means positive 
and blue negative, for that reason, among all nesting 
vectors (<#, ~ 7 r), Qi diverges first. The divergence at 
this momentum stems from nesting. 

As long as we increase J7, U c ( q) remains negative in 
a smaller region, leading to the displacement of the 
divergence to Q 2 . Why Q 2 does not change with U can 
be understood from the bare susceptibility structure, 
Xo- Xo in the entire Brillouin Zone (only <75 matters) 
can be divided in three zones: 0 < < O.I 7 r/b, 

O.lir/b < qb < O. 67 r/b and 0 . 67 t /6 < q^ < 1 r/b (and 
symmetric regions). In the first zone the susceptibility 
increases sharply due to the warping, increasing q 5 
we can connect more points of the Fermi surface; in 
the second region the system only have access to the 
Fermi sheets at one side, increasing slowly the value 
of the particle-hole susceptibility; in the third region 
connections among the two pairs of sheets gives also a 
rapid enhancement. In our case, the range of U below 
the SDW ordered phase, makes the negative U c { q) to 
be in the second region of the bare susceptibility, since 
this region is slowly q-dependent, we observe a minimal 
change of Q 2 with increasing U. The divergence of 
the charge susceptibility at this momentum is due to 
interactions and the softening can be described as a 
critical mode similar to the one found in Ref. HU In 
Fig. [ 6 ] we see a critical exponent of 

The transition from Qi to Q 2 ordering phases, is 
also shown in Fig. [5| The upper panels show the 
frequency against momentum of the charge suscep¬ 
tibility (maximum eigenvalue, which is significantly 
larger than the other three), the lower panels show 
the charge susceptibility (maximum eigenvalue) at 
zero frequency. On the left hand panels U = 0.2 
while on the right panels U = 0.4. We observe a 
change in the spectral weight of the collective mode 
from Qi to q 2 when U is increased. Moreover, while 
the weight at Qi exists at any value of V, at Q 2 the 
mode softens signaling at the proximity to the transition. 

Near the SDW region we found superconductivity in 
d x 2_ y 2 channel.This behavior is consistent with the ex¬ 
pected for a quasi-one dimensional square lattice at quar¬ 
ter fil l ing . A very recent preprint!^, proposes an order 
parameter with different sign in each band and a total 
of three node planes in the b-direction (and two more in 
the c-direction) at V = 0 . We skip the bracketed descrip¬ 
tion and project the wave (here we call it f x ) with our 
methodology. The results Fig. [3] show that both d- and 
/-channels are very close, with the f x dominating. 

As long as experiments do not show signatures of SDW 
gap opening or magnetic responsd^, we can work with 
lower U values, to avoid strong spin fluctuations. The 
Coulomb interactions are comparable with 26 . 

Near the CDW or CO regions of the phase diagram 
we found triplet superconductivity in p y channel, with 
nodes at the Fermi surface. Near the Qi CDW region, 



FIG. 7: (Color online) The background represents the charge 
fluctuations near Qi and on top of that are the CO tran¬ 
sition line (green) and superconducting transition line (solid 
black) for U = 0.2 and (U — 0.4) in right (left) panel. The 
dashed black line is the CO transition due to Qi if the other 
order were not present. The blue line is the superconducting 
transition. 


we found a narrow stripe 3 -^ of superconductivity due to 
charge fluctuations at Qi. We observe Qi is a nesting 
vector connecting all the Fermi surface with different 
phase of the order parameter. See inset Fig. [3] 

From this study, apparently we can design the interac¬ 
tions Fig. [4] (c) to be minimum in a given momentum, in 
such a way that favors superconductivity with a certain 
order parameter. Nevertheless, we need to take into 
account the bare susceptibility structure and the orbital 
distribution in real space. In the present model, the bare 
susceptibility is peaked at « 7r, and the divergence 
at Qi is favored by perpendicular Coulomb interactions 
V± and W± whereas interactions along the chains does 
not distinguish momenta in the 6 -direction. 


The Q 2 momentum is not involved in the vertex 
calculation (Eq. 10) so, we are still able to work with 
the superconducting vertex since it has not divergences. 
In that region, strong charge fluctuations still persist 
at Qi, due to nesting, and superconductivity would be 
found if another charge ordered phase were not present. 


A. Coexistence in the model 

If the order parameter of the CO is small, and assuming 
that the q 2 modulation does not open a gap at Ep, we 
consider now the possibility of coexistence with SC in 
this model, even though it may not have relevance for 
the material. 

We study the coexisting region with temperature (See 
Fig. Q. The charge fluctuations have a reentrant be¬ 
havior in RPA approaclP^, due to the fact that the bare 
susceptibility Xo(<T^) is maximum in energy (cj <C £), 
when q connects different points of the band structure 
approximately uj away from Fermi level. In that case Qi 
is a nesting vector and the maximum is close in energy, 
uj = 0.01 « T. However Q 2 exhibits its peak of reentrant 
behavior at a larger energy, and we only see the decrease 
of critical V with temperature. In Fig. [ 7 ] (a) the criti- 
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V. ELIASHBERG EQUATION WITH A 
REDUCED VERTEX 


3 

2.5 
2 

1.5 
1 

0.5 
0 

“ -200 -100 0 100 200 

n E/A 0 

FIG. 8: (Color online) Left. Gap in Matsubara frequency 
(iu ; n )at T = 0.01, U = 0.2 and V = 0.676, green line. And 
Lorentizan fit, blue line. Right. Conductance against energy 
for different superconducting order parameters. We can dis¬ 
tinguish p y wave, a is the angle of the order parameter with 
the junction. Z is the height of the tunnel barrier, in this case 
the insulating phase.® 

cal momentum changes from Qi at low temperatures to 
Q 2 - We observed that change from charge susceptibility 
and it is represented by the change of behavior of the 
CO line. The temperature makes the bare susceptibil¬ 
ity softer, lowering the value of Xo(Qi) and shifting the 
critical momentum to Q 2 . 


In the previous section we have shown that the su¬ 
perconducting vertex is dominated by the charge suscep¬ 
tibility near Qi, see (Eq. ([T()|)). It results that using 
just a few vertex momenta we reproduce the X py value. 
However, we cannot reduce easily the 4 orbital model for 
a simpler tight-binding, since near Qi the bands has a 
similar weight in the four orbitals. For that reason, we 
select the larger values of the p airin g vertex (calculated 
with the 4 orbital model) (Eq. |8|9| ) and as a result, we 
see that less than 10% of the vertex is enough to get more 
than 90% of X py value. In order to reproduce X py from 
momenta near Qi we need to multiply the value of the 
vertex by 4, otherwise we need to include momenta near 
q = 0 because it is significantly large,and connects many 
pairs of Fermi surface momenta. 

Moreover, we see that the value of the pairing vertex is 
almost independent of q y . All those simplifications, allow 
us to work with a simpler model and solve the linear 
Eliashberg equation in Matsubara frequencies, given by: 


a= 0 Z= 5 




— 1 


A ra E(k ,iu> n ) = -^ G 0 (k',iw^)r triplet (k,k , ,iw n - iu' n )G° (-k',-iu' n )T,(k', iuj'J 


k 'iu' 


(ii) 


where 


G°(kX) sp = ^^ 


<( k ) <(k) 
,(k)) 


( 12 ) 


are 4 by 4 matrices, and r triplet is also a matrix defined 
in Eq. © but with iuj n dependence coming from the 
bare susceptibility (Eq. [4] in Eq. [7]). We calculate the 
momentum dependence of the gap, by projecting on the 
p y order parameter: E(k, iuj n ) = f(iuj n ) sin(fc c c). The 
result (shown in Fig. [8| can be fitted by a Lorentzian 
plus a constant; analytic continuation or Pade approx- 
imants gives the same result, a real constant for the 
relevant frequencies. Provided the gap value is small, 
only low frequencies are relevant for experiments, as 
Normal-Insulator-Superconducting junctions—. (See 
Fig. nu) 


VI. DISCUSSION AND CONCLUSIONS 

As was previously mentioned, Lio.gMo60i7 exhibits 
signatures of Luttinger liquid behavior for a wide range 
of temperatures. Thus, it is important to discuss its re¬ 
lation with the physics described above. 


In Fig. [9j we present a schematic phase diagram 
for the model and consider its relevance for the physics 
of LiPB. Merging renormalization group estimation of 
the crossover temperature, the RPA calculations for the 
CDW and considering the fluctuation exchange in the su¬ 
perconducting vertex, we compose a schematic diagram 
Fig.U At high temperatures, the metallic phase is a LL. 
As the temperature goes down, the perpendicular hop¬ 
ping drives the system through a crossover to a Fermi 
liquid and the inter-chain Coulomb interactions through 
a thermodynamic phase transition to a CDW. Our analy¬ 
sis of the SC vertex comprises the dashed horizontal line. 
Since we are working at temperatures well below Tll , 
the use of RPA, as perturbation theory of the essentially 
free electron system, is well justified as a starting point. 
In other words, we are able to describe the superconduc¬ 
tivity as an instability of a Fermi liquid, in spite of the 
normal phase of the material being a LL. On the other 
hand, the behavior of the material as the temperature 
goes down, seems to be represented by the solid vertic al 
line. This statement is based on the spectroscopies 9 10 at 
temperatures right above T c . The density of states show 
power-law behavior very similar to the ones observed at 
much higher temperatures and similar values of a. Plac¬ 
ing the material slightly on the left of that vertical line 
would imply an interesting crossover from one NFL (the 
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FIG. 9: (Color online) Schematic phase diagram for the LiPB. 
The present study comprises the green horizontal arrow, we 
believe the temperature dependence of the real material is 
represented by the vertical orange arrow. 

LL) to another NFL (FL + strong charge fluctuations). 
At a purely qualitative level, no evidence of a Fermi edge 
developing at low temperatures has been observed and 
the experimental values of a seem to increase (instead of 
decrease). However, both alternatives rely on details of 
the model and should be quantitatively contrasted with 
the spectroscopies. 

The dashed line in Fig. [9] shows the di men sional 
crossover from Luttinger liquid to Fermi Liquid^^. The 
small value of the perpendicular hopping suggests consid¬ 
ering it as a perturbation. Based on the renormalization 
group approach, we can estimate the crossover tempera¬ 
ture to be ^ t (^) 1-Q! where a is the exponent for 
the single-particle density of states. In Fig. [lO]we show 
the estimated dimensional crossover 

for Luttinger chains coupled with Hubbard and V in¬ 
teractions. The value of a is computed using interaction 
parameters U and V following the CO border shown in 
Fig. U W is set to zero. Note that the same Coulomb 
interactions driving the charge ordering, allow for large 
values of a. Therefore, we expect Tll to be very small 
when the CDW is approached . This fact opens the pos¬ 
sibility for a direct transition from the LL to the super- 
conduting phase. It would be interesting to study this 
possibility with techniques similar to those used in Ref. 

M 


The charge ordering transition for RPA apparently oc¬ 
curs at arbitrarily large temperatures as V increases, but 
we expect the slight modifications presented in Ref. 43 
which considers how the fluctuations effects modify the 
Greens function self-consistently, evaluting also vertex 
corrections. Other details like the reentrant behavior for 
the charge ordering transition, typical of RPA calcula¬ 
tions, are unessential for the physics of the system. 

To summarize, we have studied a microscopic extended 
Hubbard model for LiPB. We have characterized the cou¬ 
plings promoting SC close to different charge ordering 
patterns. A detailed analyisis within the RPA approxi- 



FIG. 10: (Color online) Dashed. Renormalization group esti¬ 
mation of the crossover temperature for the critical V values 
found for each U. Solid. The exponent in the Luttinger den¬ 
sity of states, a. 

mation of the vertex shows triplet superconductivity with 
nodes on the Fermi surface close to those ordered phases. 
The relevance of these results is discussed in terms of the 
general experimental perspective of the material. 
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